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I Abstract 

The Monte Carlo Hamiltonian method developed recently allows to investigate ground 
state and low-lying excited states of a quantum system, using Monte Carlo algorithm with im- 
portance sampling. However, conventional MC algorithm has some difficulties when applying 
to inverse potentials. We propose to use effective potential and extrapolation method to solve 
. the problem. We present examples from the hydrogen system. 
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; 1 Introduction 



Qh. There are two standard approaches in quantum theory: Hamiltonian and Lagrangian formulations. 
^ \ A comparison is given in Tab^ The Lagrangian formulation is very suitable for applying the 
Monte Carlo (MC) method to systems with many degrees of freedom, and in the last two decades, 
it has been widely applied to lattice gauge theory 01 • In the standard Lagrangian MC method, 
^ I however, it is extremely difficult to compute the spectrum and wave function beyond the ground 



' state. 



Wave functions in conjunction with the energy spectrum contain more physical information 
than the energy spectrum alone. Although lattice QCD simulations in the Lagrangian formulation 
give good estimates of the hadron masses, one is yet far from a comprehensive understanding of 
hadrons. Let us take as example a new type of hardrons made of gluons, the so-called glueballs. 
Lattice QCD calculations 01 IHl predict the mass of the lightest glueball with quantum number 
J^'-' = O"*""*", to be 1650 it WOMeV. Experimentally, there are at least two candidates: /o(1500) 
and /j(1710). The investigation of the glueball production and decays can certainly provide 
additional important information for experimental determination of a glueball. Therefore, it is 
important to be able to compute the glueball wave function. 

*LXQ is supported by the Key Project of National Science Foundation (10235040), National and Guangdong 
Ministries of Education, and Foundation of Zhongshan University Advanced Research Center. HK is supported by 
by NSERC Canada. 
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Table 1: Comparison of the conventional methods in the standard formulations. 



Formulation 


Hamiltonian 


Lagrangian 


Approach 


Schrodinger Eq. 


Path Integral 




H\En) = En\En) 


/r>\ _ j [dx\Olx\exp(-Slx\/n) 


7-j r-j- 

Algorithm 


c — '■ '■ 

series expansion, 

variational 
Runge-Kutta ... 


MC with importance sampling 


Advantage 


Both the ground state, 


It generates the most 




and the excited states 


important configs. 




can be computed. 


for the measurements. 


Disadvantage 


Analytical methods 


It is difficult to study 




are too tedious for 


the excited states. 




many body systems; 






Runge-Kutta works 






only in 1-D. 





A natural question is whether in Lagrangian MC simulations one can construct an effective 
Hamiltonian? We have recently proposed a new approach [S] (named Monte Carlo Hamiltonian 
method) to investigate this problem. We start out from the action of the theory. We compute 
the transition matrix for all transitions using MC with importance sampling, where the states are 
taken from a basis in Hilbert space. As a result, we find an effective Hamiltonian, i.e. its spectrum 
and eigen states. A lot of models pi 171 1511^1 ITUIITTIIT^IT^ in quantum mechanics (QM) have been 
used to test the method. This method has also been applied to scalar field theories |14 | II 5 | [TB I 

Although MC with importance sampling has been very successful for QM with many local 
potentials, it is very difficult for QM with inverse potentials, due to the singularity. In this paper, 
we propose effective potential and extrapolation method to solve the problem. The hydrogen 
system is used as an example. 

The rest of the paper is organized as follows. In Sect|H and Sect. 01 the basic ideas of MC 
Hamiltonian are reviewed. In Sect. |1J we describe a method to analyze the errors of the spectrum 
and wave functions. In Sect. El a method to investigate the inverse potential is presented. The 
numerical results for the hydrogen system are given in Sect. [Bland summarized in Sect. [7| 

2 Effective Hamiltonian 

Let us discuss the construction of the effective Hamiltonian in several steps First, consider in 
quantum mechanics in 1-dim the motion of a single particle of mass m under the influence of a 
local potential V{x). Its classical action is given by 

S = Jdt (^x^ - V{x)^ . (1) 
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Given the classical action, one can determine Q.M. transition amplitudes. Like in Lagrangian 
lattice field theory, we use imaginary time in what follows. We consider the transition amplitude 
for some finite time T (tin = 0, tfi = T) and for all combinations of positions Xi,Xj. Here Xj, 
Xj run over a finite discrete set of points {xi, . . . ,xn} located on the real axis. Suppose these 
points are equidistantly distributed (spacing Ax) over some interval [—L,+L]. The transition 
amplitudes are given by the (Euclidean) path integral, 



Mij(T) = J[dx]exp[-SE[x]/h] 



(2) 



where Se = J dt (mx^/2 + y(x)) denotes the Euclidean action. Suppose for the moment that 
those transition amplitudes were known. From that one can construct an approximate, i.e. effec- 
tive Hamiltonian. Note that the matrix M(T) = [Mij{T)]NxN is a real, positive and Hermitian 
matrix. It can be factorized into a unitary matrix U and a real diagonal matrix D(T), 

M{T) = D{T) U . (3) 

Because the matrix M{T) can be expressed in terms of the full Hamiltonian H by 

Mij{T) = {xi\ expi-HT/n)\xj) , (4) 

the matices U and D have the following meaning, 



Ul =< Xi\Elff > 

eff 



Dk{T)=eM-K''T/n) , (5) 

i.e., the k — th eigenvector \E'^^-^ > of the effective Hamiltonian ii^e// can be identified with the 
k — th column of matrix and the energy eigenvalues E^-^-^ of Heff can be identified with the 
logarithm of the diagonal matrix elements of D{T). This yields an effective Hamiltonian, 

N 

Hen = Y.\K'' >K'' <K''\- (6) 
fc=i 

Note that in the above we have been mathematically a bit sloppy. The states |xi) are not Hilbert 
states. We have to replace \xi) by some "localized" Hilbert state. This can be done by introducing 
box states. We associate to each Xj some box state hi, defined by 

, . _ J 1/ VAxi if Xj < X < Xi+i , . 

^iy'^'-y else ^'^ 

where Axj = .Xj+i — x,-. When wc use an equidistant distribution of Xj, i.e. Axj = Ax, we refer 
to the basis of box states as the regular basis. 
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3 Matrix elements 



We compute the matrix element Mij{T) directly from the action via Monte Carlo with importance 
sampling. This is done by writing each transition matrix element as a ratio of two path integrals. 
This is done by splitting the action into a non-interacting part and an interacting part, 

Se = So + Sv ■ (8) 

This allows to express the transition amplitude by 



M,, (T) = (T) — , (9) 



z,0 



y,T 

z,0 



Here exp(— 5o/?i) is the weight factor and ex.p{—Sv/h) is the observable. Mf^\T) stands for 



the transition amplitude of the noninteracting system, which is (almost) known analytically. For 
details see ref.lH]. Carrying out these steps allows to construct an effective Hamiltonian, which 
has turned out to reproduce well low energy physics of a lot of quantum systems [HI El El mH 

nn m US m 113 ua nzj . 



4 Error analysis 

From errors of the matrix elements AM in Monte Carlo simulations, one can obtain the errors 
of the eigen value En and wave function of the effective Hamiltonian. Suppose H is the 
Hamiltonian of the theory, H' is a perturbation term induced by error in the simulation, then we 
have the following equation for the transition matrix element 



/ iH + H')T\ ( HT\ T f HT\ , 



which implies 



or 



/ (H + H')T\ ( HT 

AM = exp(-^-^)-exp(- — 

H' w -^Ar'^AM. (12) 

According to perturbation theory, the first order correction to the eigen value and wave function 
is 

= E^r^'/'n. (13) 
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where 



(14) 



Eq. (fT^ gives an estimate of the wave function. 

5 Inverse potential: the hydrogen atom 

The Schrodinger equation of a hydrogen is 



2m 



where 



V^*(r, e, if) + y(r)^(r, 9, if) = E^{r, 9, f), 



V{r) 



Because of the spherical symmetry, the wave function can be factorized as 

^{r,e,ip)=i;{r)Yie,ip), 
and the radial wave function satisfies 



2m 



dr \ dr 



V{r) + 



l{l + l)h^ 
2mr'^ 



'ilj{r) = Etp. 



Defining a function u{r) by 

u{r) = rip{r), 

and considering the case of angular- momentum quantum number zero (/ = 0), one obtains 

d^u 



2m dr"^ 



+ V{r)u = Eu 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



which reduces to a one dimensional QM problem. This system has the following analytic result 
for the energy 



En 



1 



2(n+l)2' 



(21) 



where n = 0, 1, ... is the principal quantum number and m = l, e = l, ^ = 1 have been chosen. 
The ground state and first excited state radial wave functions are 



V'o(?-) = 2exp(-r), 

i^i{r) = (2 - r)exp(-r/2). 



(22) 
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Our purpose is to investigate the 1-d system NOT by solving Eq. (|2U() . but by MC Hamiltonian 
method. Unfortunately, due to the divergence of V{r) at r = 0, it is not convenient for MC 
simulations. 

To solve the problem, we replace V{r) by an effective potential with cut-off parameter R 

VKir) = 5Sf^§^V(.). (23) 

exp(2r//t) + 1 

Vii{r) approaches to V{r) when — > and Vfi{r) is analytic for any value of finite -R, which will 
be used in numerical simulations. 



6 Numerical results 

The coordinates for the box state (see Eq. ((7|) are chosen to be 

ri = i*dr, i = l,2,...,N. (24) 

The simulation parameters are m = 1, e = 1, h = 1, T = 0.85, dr = 0.99 and = 100. Using the 
Metropolis algorithm and method described in Sect. |31 we obtain the matrix elements M„/„ with 
potential VR(r) for some non-zero R. Then we compute the eigenvalues and eigenvectors using 
the method described in Sect. Hand analyze their errors using the method described in Sect. 01 

Fig. n plots the ground state energy of the hydrogen system with effective potential VR{r) 
at R =0.3, 0.4, 0.5, 0.6, 0.7 and 0.8. Between initial and final r^, 600 paths are used in MC 
simulations. We use the least square method to extrapolate the data to the R = limit. The 
result agrees with the analytic one within large error bar. Fig. |2l plots the first excited state 
energy. The agreement is greatly improved if more paths (1000) are used in MC simulations. 
This is shown in Fig. 13] and Fig. ^ 

Fig. El plots the ground state wave function of the hydrogen system extrapolated to i? = 0. 
Between initial and final rj, 1000 paths are used in MC simulations. Fig. Elshows the first excited 
state wave function. They are consistent with the analytic ones. 

7 Summary 

In the preceding sections, we have extended the effective Hamiltonian method to QM with inverse 
potentials, and taken the hydrogen system as an example. The results are encouraging. We believe 
that the application of the algorithm to more complicated systems will be very interesting. It will 
allow a non-perturbative investigation of physics beyond the ground state. 
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l=0,update=600,T=0.85,N=100,dr=0.99 

-0.2 I ' ^ ' ^ ' ^ ' ^ 
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R 

Figure 1: Ground state energy of the hydrogen system with effective potential Vji{r) at R =0.3, 
0.4, 0.5, 0.6, 0.7 and 0.8. The circle at i? = is the analytic result and the line is the fitting 
curve. Between initial and final r^, 600 paths are used in MC simulations. 
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Figure 2: The same as Fig. ^ but for the first excited state energy. 
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I=0,update=1 000,T=0.85,N=1 00,dr=0.99 




Figure 3: Ground state energy of the hydrogen system with effective potential Vji{r) at R =0.3, 
0.4, 0.5, 0.6, 0.7 and 0.8. The circle at i? = is the analytic result and the line is the fitting 
curve. Between initial and final r^, 1000 paths are used in MC simulations. 
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Figure 4: The same as Fig. |31 but for the first excited state energy. 
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1=0, update=1000,T=0.85, N=100,dr=0.99 
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Figure 5: Ground state wave function of the hydrogen system extrapolated to -R = 0. The hne 
is the analytical result. Between initial and final rj, 1000 paths are used in MC simulations. 
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Figure 6: The same as Fig. |21 but for the first excited state wave function. 
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